1 
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dimensions 

Akira Niegawa 

Graduate School of Science, Osaka City University, Sumiyoshi-ku, Osaka 558-8585, 

JAPAN 

We show that the coexisting phase of quark-antiquark and diquark condensed phases is 
absent in the cold quark matter in the 2+1 dimensional extended Gross-Neveu model, which 
is in sharp contrast to the case of 3 + 1 dimensional Nambu-Jona-Lasinio model. 

§1. Introduction 

Extensive studies for last decade on quark matters have disclosed the existence 
of various phases. Among those are the quark-antiquark (qq) condensed phase, ^ 
different diquark (qq) condensed phases. Furthermore, it has been suggested 2 ) that 
the coexisting phase of qq and qq condensations also exists. (For recent reviews, see, 
e.g.,. 3 )) 

One of the important theoretical models for studying this issue is the extended 
Nambu-Jona-Lasinio (NJL) model, which is an low energy effective theory of QCD. 3 ) 
The Gross-Neveu (GN) model 4 ) proposed in 1974 is the counterpart of the NJL model 
in 1 + 1 and 2 + 1 dimensional spacetime. Study of the field theoretic models in lower 
spacetime dimension are interesting in itself as in the solid-state physics. Various 
works have been devoted to the study of the GN model. (References together with 
a brief survey of them are given in. 5 )) Recently, the phase structure of the 2 + 1 
dimensional (3D) extended GN model has been studied within the mean-field ap- 
proximation. 5 ) In, 5 ) quarks are assigned to the lowest nontrivial (2-dimensional) 
representation of the 0(2, 1) group, which we refer to as the 2d-spinor quarks. The 
model contains essentially two parameters, the qq coupling constant Gs and qq cou- 
pling constant Go- Through numerical analyses for some sets of values of (Gs, Gd), 
it has been shown that there appears qq- and gg-condensed phases but no coexisting 
phase appears, which is in sharp contrast to the case of NJL model. 

The purpose of this note is to analytically prove the absence of the coexisting 
phase for any values of Gs (> 0) and Gd- We restrict to the cold "quark matter" 
with zero temperature (T = 0). 

§2. Preliminary 

The Lagrangian of the extended GN model with 2d-spinor quarks reads 



C = q(i$+ n-yo) q + G s (qq? + G D (qr 2 X 2 q c ) (q c r 2 X 2 q) , (2.1) 

typeset using VTpT^.ds (Ver.0.9) 
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where /i is the quark chemical potential. The quark fields, q and q, are the doublets 
in the flavor space and the triplets in the color space. q c and q c are the charge- 
conjugated fields. The Pauli matrix t 2 acts on the flavor space, while the Gell-Mann 
matrix acts on the color space. Employing the mean-field approximation, we get 



C = q(i0 -a + fjrfa)q- ^ A (qT 2 \ 2 q c ) ~ -A* {q c T 2 \ 2 q) 
a 2 \A\ 2 



AGs *G D 



(2.2) 



Here 



a = -2G s (qq) , A = -2G d (?t 2 \ 2 q) , (2.3) 



where a and A are the order parameters for the qq and qq condensations, respectively. 
Computation of the thermodynamic potential at T = yields 5 ^ 

rr 2 I A 1 2 f'°° ATP 

n = ir s + ^ - / i E I 2 ^ + E *> + M + *z)*J . (2.4) 

where 



= E 2 + fj 2 + \A\ 2 ± 2 V / BV+^W 1 
£ = y^s + a 2 5 E ± = E±(i. (2.5) 



Changing the integration variable E —* £ through £ = ^ E 2 + a 2 \A\ 2 /[i 2 ± /x, the 
integrals in (|2.4p that contain may be carried out analytically. We perform 
renormalization as in, 5 ) and choose, without loss of generality, A to be real and 
positive. Then, we introduce following dimensionless quantities: 

x = alii (> 0) , y = A/n (> 1) , 
A _ 2vr / 1 3a\ g 2vr / 1 a 



27T 

I2(s,y) = -jfl( < 7,-d) j (2.6) 



where a is the renormalization scale and A = \J A 2 + fi 2 . After all this, we obtain 

n(x,y) = i {(x + l) 3 + |x - 1| 3 } + i {(x- l)|x - 1| - (x+ l) 2 } 
2 

+ 3 {(« + yf + \x- y\ 3 } + (xy - l)|x - y\ - (xy + l)(x + y) 
, / 2 i\ / 2 -a, xy + l + x + y 

+ (x -l)(y -l)ln — 1 f . (2.7) 

xy — 1 + \x — y\ 

Note that, for the normal phase, (x,y) = (0, 1), (2{x = 0,y = 1) = —1. Throughout 
in this paper, we restrict to the case < A, the case in which the quark-antiquark 
condensate can appear. 
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2.1. Quark- antiquark condensation 

Here we study the case with y = 1 (A = 0): 

fi^z) = = 1) = (3 -,4)x 2 - 1. (2.8) 

Straightforward analysis yields 

^{Ia3/27 ££|I} = * < 2 ' 9 > 

2.2. Diquark condensation 

Here we study the case with x = (a = 0): 

= ^y 3 - B(y 2 -l)-2y- (y 2 - l)L(y) - I , (2.10) 
where L(y) = ln[(y + l)/(y — 1)]. & y (y) is minimum at y = yo : 

2y -L(y ) = £, (2.11) 

n v (vo) = -^(2j/o + i) • ( 2 - 12 ) 

From (|2.1ip . Max(l, B/2) < yo an d then we have, for any B, 

Q y {y Q ) < -l = f2(x = 0, y = l). (2.13) 

Thus, the diquark condensed state is energetically favored over the "normal vacuum" . 

We now find the upper bound of f2 y (y) in the region 2 < B, which we will use 
extensively in the sequel. Since L(y) monotonically decreases with increasing y and 
its curvature is positive, 



L(y) < L(B/2) + 



dL(y) 
dy 



B 

, V ~ 2 

y=B/2 * 



Substituting this into (|2.1ip . we have 



V + ^-; + L(B/2). (2.14) 



B < 2y + ^^yo ~ - L(B/2) . (2.15) 



Solving this, we obtain 

B 2 — 4 2 B 2 — 4 
VoZ^uT + B + ^BtW/*)- (2-16) 

Here and in the following, we frequently use the inequality: 

2(1/ -1) < (y 2 -l)L(y) <2y (1 < y) . (2.17) 
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Using this in (12,16p . we obtain 



B A(B — 2) , ., , 
y,>- + ^r^ (2-18) 



Substituting ([2TT8j) into ([2T2]) . we obtain 



^ y (yo) < -§ (2 y ; 3 + lj <-— + | (=fl v ) (2<S). (2.19) 



£ 3 „ 5 



We will use f2 y , an the upper bound of Q y . Numerical computation shows that 
[fiy(yo)] B=2 = —3.515..., while Q y = —1. 

In the following sections we show that f2 > Min fi y ^ > Min ^Ac, Q y ^j for 

any value of (0 <)^4 and B, so that the quark-antiquark and diquark condensates do 
not coexist . 

§3. Region I (a; < 1 < y) 

We divide the whole region of x and y into the three regions; (0 <) x < 1 < y 
(Region I), (1 <)y < x (Region II) and 1 < x < y (Region III). In this section, we 
analyze £2(x, y) in the Region I. Analyses in the Region II and Region III will be 
made in subsequent sections. 

In the region (0 <)x < 1 < y (Region I), f2(x,y) reads 

Q(x, y) = F{y)x 2 + ^ - B(y 2 - 1) - 2y - (y 2 - l)L(y) - i , (3.1) 

F(y) = l + 2y-A + (y 2 - l)L(y) . (3.2) 

For < F(y), Q(x,y) > Q{x = 0, y) = Q y {y). Then we restrict to the region 
F(y) < 0, where 

f2(x,y) > Q{x = l,y) 

= \y Z ~ B{y 2 - 1) - A + | (= QM) . (3.3) 

f2j(y) is minimum at y = B/2 and monotonically increases with increasing (B/2 < 
)y. Since 1 < y for B < 2, Qj(y) > Qj(y = 1) = 2- A> Q x . The equality holds 
for A = 3, Qi = -1. Since f2 y (y ) < -1 (Eq. ([213]) ). &i > O y {yo) there. 
For 2<B, f2(x,y) > i?/(y = B/2). We first study 

Ayfi^y = B/2) = Qi(y = B/2) - Q y = 2B - A - 1 , (3.4) 

which is nonpositive for 

^4 + 1 

(3 <) 2B - 1 < A or B < — ^— . (3.5) 

We next study 

A x fij(y = B/2) = f2 T (y = B/2) - fi x = ^ - *L + B - A+^ (3.6) 
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in the region (13. 5p . A x ilj(y = B/2) monotonically decreases with increasing B 
(for 2 < B). Then, from ([33]), A x Qi{y = B/2) > A x fi I (y = B/2) ^ A+1 2 

(= H(A)), which is a polynomial in A. It can easily be seen that H{A) increases 
with increasing A (3 < A) and < H(A). Equality holds for A = 3 and then B = 2, 
where H(A) > f2 y (yo)- 

Note that we have not used the constraint on y that comes from F(y) < (see 
above after (|3.2p ). 

§4. Region II ((1 <) y < x) 

In this region, 

(l[x, y) = 2x 3 - Ax 2 + (y 2 - l)G(x) (4.1) 
G(x) = 2x - B + {x 2 - l)L{x) . (4.2) 

For < G(x), f2(x,y) > Q(x,y = 1) = f2 x (x), and then we restrict to the region 
G(x) < 0, where 

^{x,y) > f2(x,y = x) 

= Ax 3 -{A + B)x 2 -2x + {x 2 - l) 2 L{x) +B (= Q n (x)\ . (4.3) 

f?u(x) is minimum at x = xo : 

5xo + 2(xq — l)L(xo) = A + B , (4.4) 
fin{x ) = \ [3x 2 - {A + B)x 2 + x + B-A] . (4.5) 

(1 <) xq exists for 5 < A + B. Using ([237]) in (|P|) . we obtain 

x <{A + B + 4)/9. (4.6) 

The region A + B < 5 

^From (|4.6p . for A + B < 5, xq < 1 and then 4?/j(x) monotonically increases 
with increasing (1 <) x, so that i?//(x) > f2jj(x = 1) = 2 — ^4 > 

The region 5 < A + B 

With the help of OTTjl . the condition G(x) < (see above after flM])) yields 
S > 2x + (x 2 — l)L(x) > 4x — 2, from which we have 

x<^. (4.7) 

Since 1 < x, 2 < 

In the following, putting aside the equation (|4.4p . we regard Qh{xq) in (|4.5p as 
a function of xo and show that AQu(xq) > 0. It can be shown that, in the region 
of our interest, 1 < xo < (A + B + 4)/9, the polynomial l?//(xo) monotonically 
decreases with increasing xq. 
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We recall the condition, Eq. (|4.7|) . 

(A + B + 4)/9 < (B + 2)/4 or A< (5ff + 2)/4 : In this region, «j/(s ) > ^//(^o = 
(A + 5 + 4)/9), which is a polynomial in ^ and B. A y fl n (x a = (A + B + 4)/9) 
(= zAyi?(.A, 2?)) is a polynomial in A and 5. It can be shown that, in the region 
of our interest, this function decreases monotonically in A. Since A < (5B + 2)/4, 
A y f2(A,B) > A y Q(A = (5B + 2)/4,B), which is a polynomial in B. This polyno- 
mial increases monotonically with increasing (2 <) B and is positive. After all of 
this, we have AQu(x,y) > 0. 

(B + 2)/4 < (A + B + 4)/9 or (5ff + 2)/4 < A : In this region, f2 H (x ) > ft/j(x = 
(5 + 2)/4). The polynomial A x f2 n (x = (B + 2)/4) (= ^17(^,5)) increases mono- 
tonically with increasing ^4 (in the region of our interest). Since (5B + 2)/4 < A, 
A x fi(A, B) > A X Q(A = (5 J B+2)/4, B), which is a polynomial in B. This polynomial 
increases monotonically with increasing (2 <) B and A x il(A = (55 + 2)/4,B) > 0. 
The equality holds for S = 2 and then A = 3. At (A, B) = (3, 2), Q - Q y > 0. 

§5. Region III (1 < x < y) 

In this region, 

Q{x, y) = ^x 3 +{2y-A + (y 2 - \)L{y)) x 2 

+ ^y 3 -2y- B(y 2 - 1) - (y 2 - l)L{y) . (5.1) 

With respect to x, fl{x,y) in (15. 1H is minimum at x = xq: 

x = A-2y-(y 2 - l)L(y) . (5.2) 

The following four regions should be studied: 

Region IIIA: (|5.2[) does not have a solution in the region, 1 < xq and 1 < y. 
Region IIIB: xq < 1. 
Region IIIC: (1 <) y < x . 
region HID: 1 < xq < y. 

Region IIIA and Region IIIB 

f2(x, y) in (|5.ip monotonically increases with increasing (1 <) x. Then, f2(x, y) > 
Q{x = l,y), which is the same as (j3.3j) in the case of Region I above. Then, 
Af2(x,y) > 0. 

Region IIIC 

fi(x,y) monotonically decreases with increasing x (1 < x < y). Then, f2(x,y) > 
Q{x = y,y), which is the same as f2n(y) in (|4.3p in the case of Region II above, and 
then we briefly describe. In the present case, from (|5.2p with (|2.17p . the condition 
y < %o yields 

(l<)y = x<^, (5.3) 
5 

which is the counterpart of (|4.7p . From (|5.3p we have A < 3. We recall here (|4.6p . 
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(A + B + 4)/9 < {A + 2)/5 or B < 2(2A - l)/5: In this region, fi n (x Q ) in 

satisfies Qh{xq) > Qh(xq = (A + B + 4)/9), and the function f2n(xo = {A + 
B + 4)/9) — (= A x f2(A, B)) is a polynomial in A and B. It can be shown 
that, in the region of our interest, this function decreases monotonically in B. Since 
B < 2(2A-l)/5, A x f2(A,B) > A X Q(A,B = 2(2A-l)/5), which is polynomial in A. 
This polynomial increases monotonically with increasing (3 <) A and A X Q(A,B = 
2(2A - l)/5) > 0. The equality holds at A = 3 and then B = 2, where J? - Q y > 0. 

(A + 2)/5 < (A + B + 4)/9 or 2(2A - l)/5 < B : In this region, /l y «j/(s ) > 
Ayfljj(xo = (A+2)/5) (= Z\ y i2(^4, B)), which increases monotonically with increas- 
ing B. Since 2(2A - l)/5 < 5, 4 v ft(A,B) > ^ y ft(A,B = 2(2,4 - l)/5), which is a 
polynomial in A. This polynomial increases with increasing (3 <) A and is positive. 

5.1. Region HID 

In this region, f2(x,y) > f2(xo,y) (Eq. (|5.ip ): 

0{x Q , y) = -^ + x + ^y 3 - B(y 2 - 1) - A , (5.4) 



where xo is as in (|5.2p . 

Using (|2.17p . the condition 1 < xq < y yield 

/ A\ A + l , . 

Max M , -J < y < , (5.5) 

from which we can restrict to the region A < 3 because 1 < y. 

Ignoring the explicit form of xq, we analyze the function H(xq, y). Since f2(xo, y) 
monotonically decreases with increasing (1 <)xq, 

Q(xq, y) > [2(x = y,y) = y 3 - B(y - 1) + y - A 

= f2(A,B;y). (5.6) 

For B < 2, I? monotonically increases with increasing (1 <)y, and then f2(A, B; y) > 
Q{A,B-y = 1) = 2 -A > Q x . 

For 2 < B, f2(A, B; y) is minimum at 



y = ^ (b + 7b^ 
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J?(A, B; y ) = --(B 2 - 3) (B + a/B 2 — 3 j + y B - A . (5.7) 

For y < yo (yo < y), f2(A, B;y) monotonically decreases (increases) with in- 
creasing y. From (15, 5h . y < (A + l)/4. Solving yo(B = Bo) = (A + l)/4 for B, we 
obtain 

B {A) = * + [3(A + l) 2 + 3(A + 1) + 16] , (5.8) 

which is an increasing function of (3 <)A. Then, 

" For B < Bq(A) , y < (A + l)/4, , > 

ForB (A)<B, (A + l)/4<y , 1 J 



8 



Akira Niegawa 



Region: 2 < B < B (A) (3 < A) 

We first analyze Z\ X J7 > A 3 /27 + So; Uo(Bq)). Numerical computation 
shows that A 3 /27 + £ ; yo(-Bo)) < for A < 4.173.... Then, we analyze A y Q 
for A < 4.18: 

A y f2(A, B; y ) > A y Q{A = 4.18, B; y ) 

It can easily be seen that this polynomial is positive for A < 3 < 4.18. 
Region: B (A) < B (3 < A) 
In this case, 

Q{A,B-y) > Q{A,B-y = {A + l)/A) 

= ^±H!-^(A + 1) + 1-1[(A + 1) 2 -16]B. (5.11) 

^Prom this, we see that A x fl < for 

1 ( , . 352 
B > 91(^ + 1) - 192 + 



108 V ^4 + 5 

" IU8 [91(j4 + 1) " 192] ^= Bl ^' ( 5 - 12 ) 

Numerical computation shows that, for A < Aq = 4.396... [Aq < A], B\ < Bq 
[Bq < B\]. Then, the region of our interest is as follows: 



YoiA<Aq, B (A) < B . 
For A < A, Bi(A)<B. 



(5.13) 



Now, we analyze A y Q for B\(A) < B and Bq(A) < B, which is polynomial in 
A and B. For A + 1 < \/32, this polynomial monotonically increases with increasing 
B, while, for \/32 < A + 1, this polynomial takes minimum at B = B2(A): 



B 2 (A) = -V(A + l) 2 -32. (5.14) 



32 < A + 1 : Numerical computation shows that, 

For A <A 2 = 9.84... , B 2 < B 



For^l 2 <A, B <B 2 . (5,15) 

On the other hand, one can analytically show that B 2 < B\. After all of this, we learn 
that the region (of A and B) that should be analyzed is {Aq,Bi) < (A,B). In this 
region, A y f2(A,B) > A y Q(A, Bi), which is polynomial in A. It is straightforward 
to show that this polynomial A y Q(A, Bi) is positive for A < A. 

A + 1 < a/32 : From the above analysis, we see that, for Aq < A < \/32 — 1, 
B\ < B, and, for (3 <)A < Aq, Bq < B. In the former case, A y Q is polynomial in 
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A and is positive. In the latter case 



A y Q > A y Q{B = B ) 

> A y Q(B = B'q) 



(5.16) 



where 




+ 3 + 



4.4 + 1 



16 



(< B (A)) 



(5.17) 



which holds in the region under study. One can see that he last expression in (|5.16p . 
which is polynomial in A, is positive in the region of interest. 



We have studied the phase structure of the extended 3D Gross-Neveu model 
with 2d-spinor quarks. The following results are obtained. 

(I) There does not appear the region where the quark-antiquark and diquark con- 
densations coexist. 

(II) The diquark condensed state is always energetically favored over the normal 
vacuum (see after (|2.13p ). 

(III) From (I) and (II), we see that, for a given qq coupling constant A, there is the 
critical value, B c for the qq coupling constant B: For B < B c , gg-condensed state 
is realized, while, for B c < B, the gg-condensed state is realized. 
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